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1 Introduction 

We present a view of aspects of mathematical physics, showing how the forms 
of gauge theory, Hamiltonian mechanics and quantum mechanics arise from 
a non-commutative framework for calculus and differential geometry. 

In this paper we assume that all constructions are performed in a Lie 
algebra A. One may take A to be a specific matrix Lie algebra, or abstract Lie 
algebra. If A is taken to be an abstract Lie algebra, then it is convenient to 
use the universal enveloping algebra so that the Lie product can be expressed 
as a commutator. In making general constructions of operators satisfying 
certain relations, it is understood that we can always begin with a free algebra 
and make a quotient algebra where the relations are satisfied. 

We build a variant of calculus on A by defining derivations as commuta- 
tors (or more generally as Lie products). That is, if for a fixed N in A we 
define V : A — > A by the formula 

VF = [F, N] — FN — NF 

then V is a derivation. Note that V satisfies the formulas 

1. V(F + C7) = V(F) + V(C7) 



2. V(FC7) = V(F)G + FV(G). 



There are many motivations for replacing derivatives by commutators, or 
more generally by the derivations induced by multiplication in a Lie algebra. 
In Section 2 we give a new motivation [11] in terms of the structure of classical 
discrete calculus. The idea behind this motivation is very simple. If f(x) 
denotes (say) a function of a real variable x, and f(x) = f(x + h) for a fixed 
increment h, the we can define the discrete derivative Df by the formula 
Df = (f — f)/h, and one finds that in this classical discrete calculus the 
Leibniz rule is not satisfied. Instead one has the basic formula for the discrete 
derivative of a product: 

D(fg) = D(f)g + fD(g). 

We correct this deviation from the Leibniz rule by introducing a new non- 
commutative operator J with the property that 

fJ = Jf, 

and we define a new discrete derivative in an extended algebra by the formula 

V(/) = JD(f). 

It follows at once that 

V(fg) = JD(f)g + JfD(g) = JD(f)g + fJD(g) = V(f)g + fV(g), 
and note that 

V(/) = (Jf - Jf)/h = (fJ - Jf)/h = [f, J/h]. 

Thus in the extended algebra, discrete derivatives are represented by com- 
mutators, and naturally satisfy the Leibniz rule. This mode of translation 
shows that we can regard models based on discrete calculus as a significant 
subset of non-commutative calculus based on commutators. 

In A there are as many derivations as there are elements of the algebra, 
and these derivations behave quite wildly with respect to one another. If we 
have the abstract concept of curvature as the non-commutation of derivations, 
then A is a highly curved world indeed. Within A we shall build, in a natural 
way, a tame world of derivations that mimics the behaviour of flat coordinates 
in Euclidean space. We will then find that the description of the structure of 
A with respect to these flat coordinates contains many of the equations and 
patterns of mathematical physics. 
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Note that for any A, B, C in A we have the Jacobi Identity 



[[A,B],C] + [[C,A],B] + [[B,C],A] = 0. 

Suppose that {Vj} is a collection of derivations on A, represented re- 
spectively by {Ni} so that Vj(F) = [F,Ni] for each F in A. We define 
the curvature of the collection {Vj} to be the collection of commutators 
{R ij = [N i ,N j ]}. 

Proposition. Let the family {Vj} be given as above with V«(F) = [F,Ni\. 
then 

[V t ,V,}F=[[N l ,N J },F] 

for all F in A. Hence the curvature of {Vj} measures the deviation of the 
cocatenations of these derivations from commutativity. 

Proof. First, 

V t (V J (F)) = [[F,N J },N l ], 
which becomes via Jacobi identity 

= -[[N j ,N i \,F]-[[N i ,F],N j ] 

= [[N i ,N j ],F\ + [[F,N i \,N s ]. 

Hence 

V i (V j (F)) = [[N i ,N j ],F\ + V j (V i (F)). 

Whence 

[V l ,V 3 ]F=[[N l ,N J ],F]. 
This proves the proposition. 

In the next sections we will see how these patterns interact with concepts 
of calculus and differential geometry, and with physical models. Section 2 
shows how multivariable discrete calculus can be reformulated as a calcu- 
lus of commutators. Section 3 discusses the case of one-variable and shows 
how Brownian walks and the diffusion constant for a Brownian walk appear 
naturally from the commutator equation [X, X] = Jk where J is the shift 
operator described above, X = [X, J / r] as above with temporal increment 
r, and k is a constant. We find that for scalar X the time series executes 
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Brownian motion with diffusion constant k/2. This shows that the diffusion 
constant of a Brownian process is a structural property of that process, in- 
dependent of considerations of probability and continuum limits. Section 
3 also discusses the relationship of Brownian process and the Schrodinger 
equation. The section ends with a discussion of other solutions to the equa- 
tion [X, X] = Jk where the elements X, X', X", • • • of the time series do not 
commute with one another. 

Section 4 sets up a general format for non- commutative calculus, and 
shows that in this format (made as flat as possible), Hamilton's equations 
are naturally satisfied. This gives a new way to look at the source of this 
structure, and we discuss the relationship with the classical derivation of 
Hamilton's equations from Newtonian physics. The section goes on to discuss 
curvature in terms of commutators, as described above. Section 5 discusses 
general equations of motion in this context. We take the general dynamical 
equation in the form 

dXi/dt = Qi 

where {Gi, • • • , Gd} is a collection of elements of A. We choose to write Gi 
relative to the flat coordinates via Gi = P% — A^. This is a definition of 
and dFj dXi = [F, Pi] . The formalism of gauge theory appears naturally. In 
particular, if 

V i (F) = [F,ft], 

then we have the curvature 

[V i ,V i ]F=[i2 y , J Fl 

and 

Rij = diAj djAi + [Ai, Aj\. 

This is the well-known formula for the curvature of a gauge connection. Sec- 
tion 5 goes on to discuss how other aspects of geometry arise naturally in 
this context, including the Levi-Civita connection (which is seen as a con- 
sequence of the Jacobi identity in an appropriate non-commutative world). 
The section includes a discussion of the relationships of these structures with 
classical physics and the Poisson bracket. 



4 



Section 6 takes up the theme of the consequences of the commutator 
[Xj,Xj] = gij that we have already seen to produce the Levi-Civita connec- 
tion for the generalized metric g^. Here we carry out a sharpening of the 
work of Tanimura, deriving that 

X r = G r + F rs X s + T rst X s X\ 

where G r is the analogue of a scalar field, F rs is the analogue of a gauge field 
and r rst is the Levi-Civita connection associated with g^. This decompositon 
of the acceleration is uniquely determined by the given framework. 

Section 7 revisits the Feynman-Dyson derivation of electromagnetism 
from commutator equations, showing that most of the derivation is indepen- 
dent of any choice of commutators, but highly dependent upon the choice 
of definitions of the derivatives involved. Without any assumptions about 
initial commutator equations, but taking the right (in some sense simplest) 
definitions of the derivatives we prove a significant generalization of the result 
of Feynman-Dyson. See Theorem 7.5. 

Theorem 7.5 With the above [given in Section 7] definitions of the opera- 
tors, and taking 

V 2 = d\ + dl + dl H = X x X and E = d t X, we have 

1. X = E + X x H 

2. V«if = 

3. d t H + VxE = HxH 

4. d t E-V x H = {d 2 t - V 2 )X 

We then apply this result to produce many discrete models of the Theorem. 
These models show that, just as the commutator [X, X] = Jk describes 
Brownian motion in one dimension, a generalization of electromagnetism 
describes the interaction of triples of time series in three dimensions. 



5 



Section 8 is a discussion of the Jacobi identity. We devote part of this 
section to a proof that general Poisson brackets (not assuming Hamilton's 
equations) satisfy the Jacobi identity. This is part of the thematic struc- 
ture of this paper. We are investigating the relationship of physics and its 
variables. When the variables are commutative it is a classical matter to 
have precise locations and standard coordinates. When the "variables" are 
non- commutative one gives up the notion of location in varying degrees, and 
gets the benefit of the extra mathematical structures in non-commutative 
worlds. The Poisson bracket is singular in that it is a way to produce a 
Lie algebra structure from the algebra of derivations of a commutative al- 
gebra. Thus the Poisson bracket is a key link betweem commutative and 
non- commutative worlds. This section is intended to make our story com- 
plete and to raise the question of how this connection really comes about. 
Section 9 is a diagrammatic extension of section 8. We show how, in a dia- 
grammatic framework, the Jacobi identity can be articulated, and how it can 
arise from purely combinatorial grounds. These are hints of further discrete 
physics. Section 10 is an epilogue, discussing the themes of the paper. 

Remark. This paper is essentially self-contained, and hence it is written in 
an elementary style. While there is a large literature on non-commutative 
geometry, emanating from the idea of replacing a space by its ring of func- 
tions, this paper is not written in that tradition. Non-commutative ge- 
ometry does occur here, in the sense of geometry occuring in the context 
of non-commutative algebra. Derivations are represented by commutators. 
There are relationships between the present work and the traditional non- 
commutative geometry, but that is a subject for further exploration. In no 
way is this paper intended to be an introduction to that subject. 

The following references in relation to non- commutative calculus are use- 
ful in comparing with our approach [5, 6, 8, 24]. Much of the present work 
is the fruit of a long series of discussions with Pierre Noyes, and we will be 
preparing collaborative papers on it. I particularly thank Eddie Oshins for 
pointing out the relevance of minimal coupling. The paper [23] also works 
with minimal coupling for the Feynman-Dyson derivation. The first remark 
about the minimal coupling occurs in the original paper by Dyson [3] , in the 
context of Poisson brackets. The paper [9] is worth reading as a companion to 
Dyson. In the present paper we generalize the minimal coupling to contexts 
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including both commutators and Poisson brackets. It is the purpose of this 
paper to indicate how non-commutative calculus can be used in foundations. 
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2 Discrete Derivatives Become Commutators 

Consider a discrete deriviative Df = (f(x + A) — f(x))/A. It is easy to see 
that D does not satisfy the Leibniz rule. In fact, if 

f(x) = f(x + A), 

then 

Df=(f-f)/A 

and one calculates that 

D(fg) = D(f)g + fD(g). 

In the limit as A goes to zero, / approaches / and the Leibniz rule is satisfied. 
Now define a shift operator J that satisfies the equation 

Jf(x + A) = f(x)J 

or equivalently 

Jf = fJ- 
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Note that the existence of J is accomplished by taking the commutative 
algebra C that we started with, and extending it to the free product of C 
with an algebra generated by the symbol J, modulo the ideal generated by 
fj - Jf for all / in C. 

Setting 

V = JD, 

we have 

V(fg) = JD{f)g + JfD(g) = JD{f)g + fJD((g). 

Hence 

V(/< ? ) = V(/)(? + /V(< ? ). 
The adjusted derivative V satisfies the Leibniz rule. 

In fact, this adjusted derivative is a commutator in the algebra of functions 
C, extended by the operator J: 

V(/) = J(f-f)/A = (fJ-Jf)/A. 

Hence 

V(/) = [/,J/A]. 

Note however that 

[x, J I A] = (xJ - Jx)/A = J(x + A - x)/A = J. 

Thus V(x) = J. This underlines the fact that these derivatives now take 
values in a non-commutative algebra. Note however, that if 

a>) = x (x - A)...(x - (n - 1)A), 

then 

V(x (n) ) = JD((x {n) ) = Jnx {n - 1} . 

Hence we can proceed in calculations with power series just as in ordinary 
discrete calculus, keeping in mind powers of J that are shifted to the left. 
That is, a typical power series should be expressed in terms of the falling 
powers x^> . We would define 

ex PA (x) = £^° =0 a>7n! 
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and find that 

V(exp&(x)) = Jexp&(x). 

The price paid for having the Leibniz rule restored and the derivatives ex- 
pressed in terms of commutators is the appearance factors of J on the left in 
final expressions of functions and derivatives. 

Note that we have 

V(a;) = [x, J/ A] = J, 

and that this writes discrete calculus in terms that satisfy the Leibniz rule 
with a step of size A. It would be convenient to have an operator P such 
that [x,P] = 1. Then [f,P] would formally mimic the usual derivative with 
respect to x, and we would have 

[x n ,P] = nx n ~\ 

Of course, we can simply posit such a P, but in fact, we can redefine J so 
that 

fJ = Jf 

where 

/» = /(* + J- 1 A). 

Then 

V(ar) = [x, J/ A] = J(x + J" 1 A - x)/A = 1 

and we can take 

P = J I A. 

In this interpretation, [/, P] = JDf = V/ where 

Df(x) = (f(x + J- 1 A)-f(x))/A. 

This double readjustment of the discrete derivative allows us to transfer 
standard calculus to an algebra of commutators. For physical applications 
however, there remains a difficulty in adding in a time variable t and allowing 
that all the other elements of the algebra should be functions of time. If the 
derivative with respect to time is represented by commutation with H, then 
we cannot assume that H commutes with x. For this reason we will not 
proceed in the rest of the paper via this method of double readjustment. 
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The cost for the double readjustment is that we must have a collec- 
tion of functions in the original algebra C such that one can sensibly define 
f(x) = f(x+ J -1 A). Polynomial and power series functions have such natu- 
ral extensions. For other function algebras it will be an interesting problem 
in analysis, and algebra, to understand the structure of such extensions of 
commutative rings of functions to non- commutative rings of functions. 



3 Time, Discrete Observation, Brownian Walks 
and the Simplest Commutator 

For temporal discrete derivatives there is a very neat interpretation of the 
shift operator of the previous section. Consider a time series {X, X', X", ...} 
with commuting scalar values. Let 

DX = X = J{X'-X)/t 

where r is an elementary time step (If X denotes a times series value at time 
t, then X' denotes the value of the series at time t + r.). The shift operator 
J is defined by the equation XJ = JX' where this refers to any point in 
the time series so that X^J = JX^ n+1 ^ for any non-negative integer n. 
Moving J across a variable from left to right, corresponds to one tick of the 
clock. We already know that this discrete, non-commutative time derivative 
satisfies the Leibniz rule. 

This derivative D also fits a significant pattern of discrete observation. 
Consider the act of observing X at a given time and the act of observing (or 
obtaining) DX at a given time. Since X and X' are ingredients in computing 
(X' — X)/t, the numerical value associated with DX, it is necessary to let 
the clock tick once, Thus, if we first observe X and then obtain DX, the 
result is different (for the X measurement) if we first obtain DX, and then 
observe X. In the second case, we shall find the value X' instead of the value 
X, due to the tick of the clock. 

1. Let XX denote the sequence: observe X, then obtain X. 

2. Let XX denote the sequence: obtain X, then observe X. 
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We then see that the evaluation of these expressions in the non-commutative 
calculus parallels the observational situation: 

XX = XJ(X' -X)/t = JX'(X' - X) jr 

XX = J(X' - X)X/t. 

The numerical evaluations for two such orderings are obtained by moving 
all occurrences of J all the way to the left. Thus we could write 

\J m A\ = A 

for an expression where A has no appearance of J. Then 

\XX\ =X'{X' -X)/t 

and 

\XX\ = (X' - X)X/t. 

Elsewhere [11] we have called this interpretation of the temporal discrete 
derivative the "Discrete Ordered Calculus" or DOC for short. 

The commutator [X, X\ expresses the difference between these two orders 
of discrete measurement. In the simplest case, where the elements of the time 
series are commuting scalars, we have 

[X, X] = XX - XX = XJ(X' -X)/t- J((X' - X) /t)X 

= J[X\X' - X) - (X' - X)X\/t = J(X' - Xfjr. 
Thus we can interpret the equation 

[X, X] = Jk 

(k a constant scalar) as 

{X 1 - Xfjr = k. 
This means that the process is a Brownian walk with spatial step 



A = ±Vkr 
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where A; is a constant. In other words, we have 

k = A 2 /r. 

We have shown that a Brownian walk with spatial step size A and time step 
r will satisfy the commutator equation above exactly when the square of the 
spatial step divided by the time step remains constant. This means that a 
given commutator equation can be satisfied by walks with arbitrarily small 
spatial step and time step, just so long as these steps are in this fixed ratio. 

Remarkably, we can identify the constant k/2 as the diffusion constant 
for the Brownian process. To make this comparison, lets recall how the 
diffusion equation usually arises in discussing Brownian motion. We are 
given a Brownian process where 

x(t + r) = x(t) ± A 

so that the time step is r and the space step is of absolute value A. We regard 
the probability of left or right steps as equal, so that if P(x, t) denotes the 
probability that the Brownian particle is at point x at time t then 

P(x, t + r) = P(x- A, t) /2 + P(x + A, t) /2. 

From this equation for the probability we can write a difference equation for 
the partial derivative of the probability with respect to time: 

(P(x, t + r)- P(x, t))/r = (h 2 /2r) [(P(x - A, t) - 2P(x, t) + P(x + A))/A 2 ] 

The expression in brackets on the right hand side is a discrete approximation 
to the second partial of P(x, t) with respect to x. Thus if the ratio C = A 2 /2r 
remains constant as the space and time intervals approach zero, then this 
equation goes in the limit to the diffusion equation 

dP(x,t)/dt = Cd 2 P(x,t)/dx 2 . 

C is called the diffusion constant for the Brownian process. 

The appearance of the diffusion constant from the observational commu- 
tator shows that this ratio is fundamental to the structure of the Brownian 
process itself, and not just to the probabilistic analysis of that process. 
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3.1 Planck's Numbers, Schrodinger's Equation and the 
Diffusion Equation 

First recall the Planck Numbers, ft is Planck's constant divided by 2tt. c is 
the speed of light. G is Newton's gravitational constant. The Planck length 
will be denoted by L, the Planck time by T and the Planck mass by M. Their 
formulas are 

M = ^ftc/G 
L = ft/Mc 
T = h/Mc 2 . 

These portions of mass, length and time are constructed from the values of 
fundamental physical constants. They have roles in physics that point to 
deeper reasons than the formal for introducing them. Here we shall see how 
they are related to the Schrodinger equation. 

Recall that Schrodinger's equation can be regarded as the diffusion equa- 
tion with an imaginary diffusion constant. Recall how this works. The 
Schrodinger equation is 

ihdtp/dt = Hi/; 

where the Hamiltonian H is given by the equation H = p 2 /2m + V where 
V(x, t) is the potential energy and p = {ft/i)d/dx is the momentum operator. 
With this we have p 2 /2m = {-ft 2 /2m)d 2 / dx 2 . Thus with V(x,t) = 0, the 
equation becomes iftd^/dt = (—ft 2 /2m)d 2 il)/dx 2 which simplifies to 

Oip/dt = (ift/2m)d 2 ^/dx 2 . 

Thus we have arrived at the form of the diffusion equation with an imagi- 
nary constant, and it is possible to make the identification with the diffusion 
equation by setting 

ft/m = A 2 /t 

where A denotes a space interval, and r denotes a time interval as explained 
in the last section about the Brownian walk. With this we can ask what 
space interval and time interval will satisfy this relationship? Remarkably, 
the answer is that this equation is satisfied when m is the Planck mass, A is 
the Planck length and r is the Planck time. For note that 

L 2 /T = (ft/Mc) 2 /(ft/Mc 2 ) = ft/M. 
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What does all this say about the nature of the Schrodinger equation it- 
self? Interpreting it as a diffusion equation with imaginary constant suggests 
comparing with the DOC equation 

[X,X] = iJk 

for a real constant k. This equation implicates a Brownian process where 
X' = X ± T where T 2 /t = ik. We can take T = \JiA where A is a real step- 
length. This gives a Brownian walk in the complex plane with the correct 
DOC diffusion constant. However, the relationship of this walk with the 
Schrodinger equation is less clear because the ip in that equation is not the 
probability for the Brownian process. To see a closer relationship we will 
take a different tack. 

Consider a discrete function ip(x,t) defined (recursively) by the following 
equation 

il>(x, t + r) = {i/2)^{x - A, t) + (1 - i)i[>(x, t) + {i/2)^{x + A, t) 

In other words, we are thinking here of a random "quantum walk" where the 
amplitude for stepping right or stepping left is proportional to i while the 
amplitude for not moving at all is proportional to (1 — i). It is then easy to 
see that ip is a discretization of 

d^/dt = {iA 2 /2r)d 2 i)/dx 2 . 

Just note that tp satisfies the difference equation 

(i/>(x, t + t) -i>(x, t))/r = (iA 2 /2T)^(x — A,t)— 2ip(x, t) +^(x + A, t))/A 2 

This gives a direct interpretation of the solution to the Schrodinger equation 
as a limit of a sum over generalized Brownian paths with complex amplitudes. 
We can then reinterpret this in DOC terms by the equation [X, X] = J(A 2 /r) 
or [X, X] = 0, each of these contingincies happening probabilistically. For a 
different (and deeper) relationship between Brownian motion and quantum 
mechanics see [25]. 
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3.2 DOC Chaos 



Along with the simple Brownian motion solution to the one dimensional com- 
mutator equation, there is a heirarchy of time series that solve this equation, 
with periodic and chaotic behaviour. These solutions can be obtained by 
taking 

X = J n Y 

where Y is a numerical scalar, and taking the commutator equation to be 

[X,X] = J 2n+1 k 

where A; is a scalar. Expanding this equation, we find 

XJ{X' - X) - J(X' - X)X = J 2n+1 k 

J n YJ(J n Y' - J n Y) - J(J n Y' - J n Y)J n Y = J 2n+1 k 

j2n+iY n + 1 (Y l _ y) — j 2n + 1 (Y n + 1 _ Y n )Y = J 2n+l k 

yn+i(y' -Y)- {Y n+l - Y n )Y = k 

Y n+l (Y' -2Y) = k-Y n Y 

Y n+1 — (k — Y n Y)/(Y' - 2Y). 

This last equation expresses the time series recursively where Y refers to the 
value of the series that is n time steps back from Y n . The first case of this 
recursion is 

Y" = (k-Y'Y)/(Y' -2Y). 

Next case is 

Y"' = {k-Y"Y)/{Y' -2Y). 

These recursions depend critically on the value of the parameter k. In the 
first case one sees periodic oscillations that (for appropriate values of k) 
destabilize and blow up, alternating between an unbounded phase and a 
bounded semi-periodic phase. In Figures 1 and 2 we illustrate the case of 
the equation 

Y" = (k-Y'Y)f(Y' -2Y) 

for k = .0001 in Figure 1 (a bounded phase) and k = .009 in Figure 2 (an 
unbounded phase). There is intricate recursive structure in this hierarchy 
and it deserves further study. 
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Figure 1 - Y~" = (.0001 - Y'Y)/{Y' - 2Y) 




Figure 2 — Y" — (.009 - Y~'Y~)/(Y~' - 2Y) 
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4 Non Commutative Calculus and Hamilton's 
Equations 



We now set up a framework for non- commutative calculus in an arbitrary 
number of dimensions. We shall assume that each derivative is represented 
by a commutator, and that the basic space and time derivatives commute 
with one another as is customary for the flat space of standard multi- variable 
calculus. This production of a flat space for calculus forms a reference domain 
within the containing Lie algebra A. 

Since all derivatives are represented by commutators, this includes the 
time derivative as well. We shall assume that there is an element H in A 
representing the time derivative. This means that 

dA/dt = [A, H] 

for any A in A. Note that it follows at once from this choice that H itself is 
time independent, since dH/dt = [H, H } = 0. We shall see that H behaves 
formally like the Hamiltonian operator in classical mechanics. 

We will assume that there is a set of coordinates {X ± , ...,X d } that are as 
flat as possible. It is assumed that the X^ all commute with one another, and 
that the derivatives with respect to them commute with one another. The 
partial derivatives with respect to Xi will be represented by a set of elements 
{P ir --,P d } with 

d t F = dF/dX l = [F,P l ] 

for any F in A. 

Since we want the equation 

dXi/dXj = S i:j , 

we need the equation 

[Xi, Pj] = 5ij. 



17 



Since we want 

didj = djdi, 

and since (as we compute in the introduction and in the next section) 

didjF - djdiF = [di, dj]F = [[P^ Pj], F], 

we see that these partial derivatives will commute with one another exactly 
when [Pj, Pj] belongs to the center of the algebra A for all choices of % and j. 

For simplicity, we shall assume that 

[P,P,]=0. 

With these choices the flat coordinates satisfy: 

[X i ,X j ] = 

[Pi,Pj}=0 
[Xi, Pj] = 5ij. 

Note that we also have 

d t F = dF/dP l = [X h F] 

so that 

d j X i = dX i /dP j = [X j ,X i \=Q 

and 

d j P i = dP i /dP j = [X j ,P^=8 ij . 

This formalism looks like bare quantum mechanics and can be so inter- 
preted, (if we take ihdA/dt = [A, H] and H the Hamiltonian operator). But 
these coordinates can also be viewed as the simplest flat set of coordinates 
for referring the description of temporal phenomena in a non-commutative 
world. There are various things to note. For example 

dPi/dt = [P, H] = -[H, P] = -dH/dX l 
dXi/dt = [X i} H] = dH/dPi. 
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Thus 



dPi/dt 



dH/dXi 



dXi/dt = dH/dPi. 



These are exactly Hamilton's equations of motion. The pattern of Hamilton's 
equations is built into the system! 



4.1 Hamilton's Equations in Classical Mechanics 

It is worth recalling how Hamilton's equations appear in classical mechanics. 
For simplicity, we shall restrict to one spatial variable q (the analog of the 
operator X) and one momentum variable p (the analog of the operator P). 
In classical mechanics in one space and one time dimension, we have the 
equations 



where the first equation is the definition of momentum of a particle of mass 
to, the second equation is the expression for the energy of the system as the 
sum of the kinetic energy p 2 /2m and the potential energy V(q). The third 
equation is Newton's law of motion. 

We see that 



p = mdq/dt, 
H = p 2 /2m + V(q), 
md 2 q/dt 2 = -dV/dq, 



dH/dp 



p/m = dq/dt 



and 



dH/dq = dV/dq 



md 2 q/dt 2 = —dp/dt. 



Thus 



dH/ dp = q 



and 



dH/dq = —p. 



These are Hamilton's equations of motion. 
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Hamilton went on to observe that for any function F of q and p, 
F = dF/dt = dF/dq q + dF/dpp 
= dF/dq OH I dp - OF /dp dH/dq. 

Thus 

F = {F,H}, 

where {F, H} is the Poisson Bracket defined by the equation 

{A, B} = dA/dqdB/dp - OA/dp dB/dq. 

Remarkably, the Poisson Bracket satisfies the Jacobi identity, and hence gives 
a Lie algebra structure on the commutative space of functions of position and 
momentum. We have shown in this section that the pattern of Hamilton's 
equations is inherent in the Lie algebra context. We shall have more to say 
about Poisson Brackets and Hamilton's equations in Section 5. 



4.2 Curvature 

Note that for any A, B, C in A we have the Jacobi Identity 

[[A,B],C\ + [[C,A],B] + [[B,C\,A]=0. 

Suppose that {Vj} is a collection of derivations on A, represented re- 
spectively by {Ni} so that Vi(F) = [F, iVj] for each F in A. We define 
the curvature of the collection {Vj} to be the collection of commutators 
{R ii = [N i ,N j \}. 

Proposition 4.1 Let the family {Vj} be given as above with Vj(-F) = 
[F, . then 

\7i,V j ]F=[[N i ,N j ],F] 

for all F in A. Hence the curvature of { Vj} measures the deviation of the 
cocatenations of these derivations from commutativity. 
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Proof. First, 

V l (V J (F)) = [[F,N J ],N t ], 
which becomes via Jacobi identity 

= -[[N J ,N l ],F]-[[N u F],N J ] 

= [[N i ,N s ],F\ + [[F,N i \,N j ]. 

Hence 

V,(V,(F)) = [[A^],F] + V^F)). 

Whence 

[V t ,V 3 ]F=[[N l ,N J ],F]. 
This proves the proposition. □ 

This elementary notion of curvature for a collection of derivations just 
measures the extent to which they do not commute with one another. We 
have a collection X of elements X of the algebra, and corresponding deriva- 
tions V x where V xZ = Z, P x for a corresponding set of elements P x repre- 
senting these derivations. This is exactly the situation of the main framework 
of this section where 

[X t ,X 3 ]=0 

[P»Pi\=0 
[Xi, Pj] = Sij 

so that 

P Xi = Pi- 

In this case we have ViZ = [Z,Pi], and the curvature of the collection is 
the collection of commutators {Rij = [Pi,Pj]}. More generally, the elements 
of the collection X may not commute with one another. In this case we 
shall define the curvature as an operator R(X, Y) defined on X x X by the 
equation 

R(X, Y)Z = [Vx, V Y ]Z - V [X , Y] Z 

in direct analogy with the usual definition in standard differential geometry. 
However, in order to do this we shall need a collection X closed under the 
commutator, so that V[x,y] is defined. The next paragraphs outline one way 
to accomplish this end. 
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We will build the notion of curvature in terms of a general concept of 
covariant differentiation. Let X' denote a specific collection of elements of 
the algebra A that we shall refer to as the variables in A. We include in A a 
notion of time in the sense that there is the temporal derivative F = [F, H] 
and we designate a single special variable T to correspond to this temporal 
derivative. In general, to each variable X there is associated a derivative Vx 
with associated action 

V X (F) = [F,P X ] 

where Px represents this derivative. (Thus we take Pt — H.) In this gener- 
ality, we make no assumptions about the commutativity of the variables or 
of the corresponding elements that represent the derivatives. 
We define the following set S of "scalar functions over X n . 

S = Center{f G A\ [X, /] = 0, X e X'}. 

That is, S is the set of / that commute with each other and with all the 
elements of X'. We then consider the closure of X' under addition and mul- 
tiplication by elements of S. Call this closure X . Define Px+fY = Px + fPy 
for X, Y e X, f e S. 

For / e S and X E X, define X[f] by the formula 

X[f] = V x (f) = [f,Px}. 

Note that even though / commutes with all the elements of X, it can still 
have non-trivial derivations with respect to these variables. 

Lemma 4.2 We have the following properties: 

1. For all X, Y e X, V x +y = V x + Vy. 

2. Forall/e<S,X,y e X,W fx (Y) =fV x (Y). 

3. For all feS,X,YeX, V x (fY) = X[f]Y + fV x (Y). 
Proof. For the first property, note that 

V x+ y{Z) = [Z,P x+ y] = [Z,P x + Py}=Vx(Z) + Vy(Z). 
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For the second property, we have 

V fx (Y) = [Y,P fX ] = [YJP X ] = [Y,f]P x +f[Y,P x ] = 0+fV x (Y) = fV x (Y). 

For the third property, note that 

V x (fY) = [fY, P x ] = [/, P X ]Y + f[Y, P x ] = X[f]Y + fV x (Y). 

This completes the proof. □ 

Definition. We define the curvature as a function R : X x X x X — > A by 
the formula 

R(X, Y)Z = [Vx, Vy]Z - V [XX] Z. 

Thus for given elements X,Y G X, the curvature operator R(X, Y) measures 
the non-commutativity of the operators V x and Vy in relation to the non- 
commutativity of X and Y. If X and Y commute, then 

R(X,Y)Z= [V X ,V Y ]Z, 

and we are returned to our initial definition of curvature for a collection of 
derivations. 



5 General Equations of Motion 

Given a set of coordinates {X i: X 2 , ■ ■ ■ X d } and dual coordinates {Pi, P 2 , • • • Pd} 
as in the previous section, a general description of dXi/dt takes the form of 
a system of equations 

dXjdt = Q i 

where {Qi, • • • , Qd} is a collection of elements of A. If we choose to write Qi 
relative to the flat coordinates via Qi = Pi — A\ (this is a definition of Aj) 
then the formalism of gauge theory appears naturally. For example, if 

Vi(F) = [F,Qi\, 

then we have the curvature 

[V J ,V,]F=[ J R lJ ,F] 
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where 

Rij = [Gi, Gj] 
= [Pi — Mi Pj — Aj] 
= -[P i ,A j ]-[A i ,P j \ + [A i ,A j ] 

= diAj djAi + [Ai, Aj\. 

This is the well-known formula that expresses the gauge field as the curvature 
of the gauge connection. From this point of view everything comes naturally 
from the assumption that all derivatives are represented by commutators, 
and that one refers all equations to the flat background coordinates. 

5.1 Curvature and Connection at the Next Level 

The dynamical law is 

dX t /dt = Xi = Pi - Ai = Gi- 
This gives rise to new commutation relations 

[X t , XA = [X h P 3 ] - [X t , A,} = Sij - dAj/dPt = 9ij 
where this equation defines g^-, and 

[Xt, Xj] = = - djAi + [Ai, Aj\. 
We define the "covariant derivative" 

V,F = [F, P - A,] = di(F) - [F, A,] = [F, X t ], 
while we can still write 

d t F = [X t ,F}. 
24 



It is natural to think that gij is analogous to a metric. This analogy is 
strongest if we assume that 

[Xi,g jk ] = 0. 

By assuming that the spatial coordinates commmute with the metric coeffi- 
cients we have that 

[Xi,g jk ] + [X h g jk \ = 0. 

Hence 

Here, we shall let 

not assuming this commutator to vanish. Then 

Ciji = [Xi,gjk] + [Xi,gjk] = digjt — Vigjk- 

A stream of consequences then follows by differentiating both sides of the 
equation 

9ij = [X h Xj\. 

We will detail these consequences in section 6. For now, we show how the 
form of the Levi-Civita connection appears naturally in this context. 
In the following we shall use D as an abbreviation for d/dt. 

The Levi-Civita connection 

T ijk = (1/2)(V ig jk + Vjgik - VkQij) 

associated with the g^ comes up almost at once from the differentiation pro- 
cess described above. To see how this happens, view the following calculation 
where 

d i d j F=[X i ,[X j ,F\]. 
We apply the operator fydj to the second time derivative of X k . 
Lemma 5.1 Let T ijk = (l/2)(V^ fc + Vjg ik - VkQij)- Then 

T ijk = (l/2)(d i d j X k - dig jk ), 
where c ijk = [X i: g jk ] = dig jk . 
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Proof. Note that by the Leibniz rule 

D([A,B]) = [A,B] + [A,B], 

we have 

9jk — [Xj,X k ] + [Xj,X k ]. 

Therefore 

didjX k = [Xi, [Xj,X k ]] 
= [Xi,g'jk - [Xj,X k ]\ 
— [Xi-,gjk] — [Xi, [Xj,X k ]} 
— [Xi,gjk\ + [X k , [Xi,Xj\] + [Xj, [X k ,Xi] 
= Ciji — [Xi, gj k ] + [X k , [Xi, Xj]] + [Xj, [X k , Xi] 

= Cij k + ViQjk — ^kdij + Vj^ifc 

Cij k 2r&jj. 

This completes the proof. □ 

It is remarkable that the form of the Levi-Civita connection comes up directly 
from this non-commutative calculus without any apriori geometric interpre- 
tation. 

The upshot of this derivation is that it confirms our interpretation of 

9i j = [X h X 3 ] = [X,, Pj\ - [X,, Aj\ = 5 tJ - dAj/dPi 

as an abstract form of metric (in the absence of any actual notion of dis- 
tance in the non-commutative world). This calls for a re-evaluation and 
reconstruction of differential geometry based on non-commutativity and the 
Jacobi identity. This is differential geometry where the fundamental concept 
is no longer parallel translation, but rather a non-commutative version of a 
physical trajectory. This approach will be the subject of a separate paper. 
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At this stage we face the mystery of the appearance of the Levi-Civita 
connection. There is a way to see that the appearance of this connection is 
not an accident, but rather quite natural. We shall explain this point of view 
in the next subsection where we discuss Poisson brackets and the connection 
of this formalism with classical physics. On the other hand, we have seen 
in this section that it is quite natural for curvature in the form of the non- 
commutativity of derivations to appear at the outset in a non-commutative 
formalism. We have also seen that this curvature and connection can be 
understood as a measurement of the deviation of the theory from the "flat" 
commutation relations of ordinary quantum mechanics. Electromagnetism 
and Yang-Mills theory can be seen as the theory of the curvature introduced 
by such a deviation. On the other hand, from the point of view of metric 
differential geometry, the Levi-Civita connection is the unique connection 
that preserves the inner product defined by the metric under the parallel 
translation defined by the connection. We would like to see that the formal 
Levi-Civita connection produced here has this property as well. 

To this end lets recall the formalism of parallel translation. The infinites- 
imal parallel translate of A is denoted by A' = A + 5 A where 

6A k = _ r k A i dX j 

where here we are writing in the usual language of vectors and differentials 
with the Einstein summation convention for repeated indices. We assume 
that the Christoffel symbols satisfy the symmetry condition rf- = Y k - V The 
inner product is given by the formula 

< A,B >= Qlj A l B 3 

Note that here the bare symbols denote vectors whose coordinates may be 
indicated by indices. The requirement that this inner product be invariant 
under parallel displacement is the requirement that 8{gijA l A>) = 0. Calcu- 
lating, one finds 

SfaA'Al) = {d k9ij )A i A^dX k + g l3 5(A*)Ai + g t ,A%A^) 
= (d k g tj )A*AJdX k - g%j r rs A r dX s Ai - g ir \ : l^Vd.X* 
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= (d k g ij )A i A^dX k - 9ij r rs A r A^dX s - g^A 1 A r dX s 
= (d k g t ,)A l ^dX k - g^-AW'd.X 1 - gJ^.AW'd.X 1 

Hence 

(d k gij) = g S j^tk + 9is^ s jk - 

From this it follows that 

^ijk = 9is^ s jk = (V 2 )( , 9fcfi'ij - dig jk + dj(g ik )). 

Certainly these notions of variation can be imported into our abstract con- 
text. The question remains how to interpret the new connection that arises. 
We now have a new covariant derivative in the form 

V l X ] = diX j + T J kl X k . 

The question is how the curvature of this connection interfaces with the 
gauge potentials that gave rise to the metric in the first place. The theme 
of this investigation has the flavor of gravity theories with a qauge theoretic 
background. We will investigate these relationships in a sequel to this paper. 

5.2 Poisson Brackets and Commutator Brackets 

Dirac [7] introduced a fundamental relationship between quantum mechanics 
and classical mechanics that is summarized by the maxim replace Poisson 
brackets by commutator brackets. Recall that the Poisson bracket {A, B} is 
defined by the formula 

{A, B} = (dA/dq)(dB/d P ) - (dA/dp)(dB/dq), 
where q and p denote classical position and momentum variables respectively. 

In our version of discrete physics the noncommuting variables are func- 
tions of time, with the time derivative itself a commutator. With 

DF = [F,J/r], 
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it follows that 

D([A,B}) = [DA, B] + [A, DB] 

for any expressions A, B in our Lie Algebra. A corresponding Leibniz rule 
for Poisson brackets would read 

(d/dt){A, B} = {dA/dt, B} + {A, dB/dt}. 

However, here there is an easily verified exact formula: 

(d/dt){A, B} = {dA/dt, B} + {A, dB/dt} - {A, B}(dq/dq + dp/ dp). 

This means that the Leibniz formula will hold for the Poisson bracket exactly 
when 

(dq/dq + dp /dp) = 0. 

This is an integrability condition that will be satisfied if p and q satisfy 
Hamilton's equations 

q = dH/ dp, 
p = —dH/dq. 

This means that q and p are following a principle of least action with re- 
spect to the Hamiltonian H. Thus we can interpret the fact D([A,B]) = 
[DA, B] + [A, DB] in the non- commutative context as an analog of the prin- 
ciple of least action. Taking the non-commutative context as fundamental, 
we say that Hamilton's equations are motivated by the presence of the Leib- 
niz rule for the discrete derivative of a commutator. The classical laws are 
obtained by following Dirac's maxim in the opposite direction! Classical 
physics is produced by following the correspondence principle upwards from 
the discrete. 

In making this backwards journey to classical physics we see how our 
earlier assertion that bare mechanics of commutators can be regarded as the 
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background for the coupling with other fields (as in the description of formal 
gauge theory), fits with Poisson brackets. The bare Poisson brackets satisfy 

{q h qj} = 

{PhPj} = o 
{QhPj} = Sij. 

In our previous formalism, we would identify Xi as the correspondent with 
qi and Pj as the correspondent of pj. And, given a classical vector potential 
A, we could write the coupling dqi/dt = pi — Ai to describe the motion 
of a particle in the presence of an electromagnetic field. Similar remarks 
apply to the analogs for gauge theory and curvature. In particular it is of 
interest to see that our derivation of the Levi-Civita connection corresponds 
to the motion of a particle in generalized coordinates that satisfies Hamilton's 
equations. The fact that such a particle moves in a geodesic according to the 
Levi-Civita connection is a classical fact. Our derivation of the Levi-Civita 
connection, interpreted in Poisson brackets, reproduces this result. 

To see how this works, let ds 2 = g^dxidxj denote the metric in the 
generalized coordinates x k - Then the velocity of the particle has square v 2 = 
(ds/dt) 2 = g^XiXj. The Lagrangian for the system is the kinetic energy 
L = mv 2 /2 = mg^XiXj/2. Then the canonical momentum is pj = dL/dxj, 
and with = Xi we have the Poisson brakets 

Sij = tii,Pj} = {xi,dL/dxj} = {xi,mg 3k x k }. 

Taking m — 1 for simplicity, we can rewrite this bracket as 

{xi, x j} = gij. 

This, in Poisson brackets, is our generalized equation of motion. 

The classical derivation applies Lagrange's equation of motion to the sys- 
tem. Lagrange's equation reads 

d/dt(dL/dxi) = dL/dxi. 

Since this equation is equivalent to Hamilton's equation of motion, it follows 
that the Poisson brackets satisfy the Leibniz rule. With this, we can proceed 
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with our derivation of the Levi-Civita connection in relation to the accelera- 
tion of the particle. In the classical derivation, one writes out the Lagrange 
equation and solves for the acceleration. The advantage of using only the 
Poisson brackets is that it shows the relationship of the connection with the 
Jacobi identity and the Leibniz rule. 

This discussion raises further questions about the nature of the general- 
ization that we have made. Originally Hermann Weyl [30] generalized clas- 
sical differential geometry and discovered gauge theory by allowing changes 
of length as well as changes of angle to appear in the holonomy. Here we 
arrive at a similar situation via the properties of a non-commutative discrete 
calculus of observations. 



6 Consequences of the Metric 

In this section we shall follow the formalism of the metric commutator equa- 
tion 

[X, h Xj] = Qij 

very far in a semi-classical context. That is, we shall set up a non-commutative 
world, and we shall make assumptions about the non-commutativity that 
bring the operators into close analogy with variables in standard calculus. In 
particular we shall regard an element F of the Lie algebra to be a "function 
of the Xi " if F commutes with the X^ and we shall assume that if F and 
G commute with the JQ, then F and G commute with each other. We call 
this the principle of commutativity. With these background assumptions, it 
is possible to get a very sharp result about the behaviour of the theory. In 
particular, the results of this section sharpen the work in [28], where special 
orderings and averages of orderings of the operators were needed to obtain 
analogous results. 

We assume that 

[Xi,Xj] = gij 
[X i ,X j ] = 
[Xi,g jk ] = 
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[9ij,9ki] = 0. 
We assume that there exists a with 

9 lj 9jk = $1 = 9i 3 9 jk = Si 

We also assume that if 

[A,Xi] = 

and 

[B,Xi]=0 

for all i, then 

[A,B]=0 

for all expressions A and B in the algebra under consideration. To say that 
[A, Xj\ = is to say the analogue of the statement that A is a function only 
of the variables Xj and not a function of the Xj. This is a stong assumption 
about the algebraic structure, and it will not be taken when we look at strictly 
discrete models. It is, however, exactly the assumption that brings the non- 
commutative algebra closest to the classical case of functions of positions and 
momenta. 

The main result of this section will be a proof that 

X r = G r + F rs X s + T rst X s X\ 

and that this decompositon of the acceleration is uniquely determined by the 
given framework. Since 

F" = [X>;X']=g ri g'tF ij , 

we can regard this result as a description of the motion of the non- commutative 
particle influenced by a scalar field G r , a qauge field F rs , and geodesic mo- 
tion with respect to the Levi-Civita connection corresponding to gij. Let us 
begin. 

Note that, as before, we have that = g^ by taking the time derivative of 
the equation [Xj,Xj] = 0. 

Note also that the Einstein summation convention (summing over repeated 
indices) is in effect when we write equations, unless otherwise specified. 
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As before, we define 

diF=[F,Xi] 

and 

diF=[Xi,F]. 

We also make the definitions 
and 

d l F = [F,X*]. 

Note that we do not assume the existence of a variable X 3 whose time deriva- 
tive is X 3 . Note that we have 

X k = 9kiX % . 
Note that it follows at once that 

diQjk = digj k 

by differentiating the equation [Xi,gj k ] = 0. 

We assume the following postulate about the time derivative of an element 
F with [Xt, F] = for all k : 

This is in accord with the concept that F is a function of the variables X«. 
Note that in one interpretation of this formalism, one of the variables Xi 
could be itself a time variable. In the next section, we shall return to three 
dimensions of space and one dimension of time, with a separate notation for 
the time variable. Here there is no restriction on the number of independent 
variables X^. 

We have the following Lemma. 

Lemma 6.1 

1. [X h Xi] = 5i. 

2. d r (g i3 )g jk + g i3 d r (g jk ) = 0. 

3. [X r , dig jk ] = 0. 
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Proof. 

[Xi,Xi\ = [X t ,g^ k X k ] = [X h g^\X k + g^[X h X k \ 
= g jk [X h X k ] = g ]k g ik = g ]k g kl = Sf. 
The second part of the proposition is an application of the Leibniz rule: 

= d r (5i) = d r (g ij g jk ) = d r (g ij )g jk + g ij d r (g jk ). 

Finally, 

[X r , digj k ] = [X r , \gjk,Xi]] = —[Xi, [X r ,gj k ]] — [g jk , [Xi,X r ]] 

= -{X i ,0] + [g jk ,g tr } =0 + = 0. 
This completes the proof of the Lemma. □ 

It follows from this lemma that d l can be regarded as d/dXi. 

We have seen that it is natural to consider the commutator of the veloci- 
ties Rij = [Xi, Xj] as a field or curvature. For the present analysis, we would 
prefer the field to commute with all the variables X k in order to identify it as 
a "function of the variables X k " . We shall find, by a computation, that -Re- 
does not so commute, but that a compensating factor arises naturally. The 
result is as follows. 

Proposition 6.2 Let F rs = [X r , X s ] + (d r g ks - d s g kr )X k and F rs = [X r , X s }. 
Then 

1. F rs and F rs commute with the variables X k . 

2. F rs = g ri g sj F i:j . 

Proof of Proposition. 

We begin by computing the commutator of X; L and R rs = [X r ,X s ] by 
using the Jacobi identity. 

[Xi, [X r ,X s ]] = —[X s , [Xi,X r ]] — [X r , [X s ,Xi\] = d s g ir — d r g is . 

Note also that 

[Xi,d r g ks ] = [Xi, [g ks ,X r ]] = —[X r , [Xi,g ks \] — [g ks , [X r ,Xi\] 
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= -[X r , [Xi,g ks ]] + [gks,9ir] = 0. 

Hence 

[Xi, (d r g ks — d s g kr )X k ] = d r g is — d s g ir . 

Therefore 

[Xi, F rs ] = [Xi, [X r , X s ] + (d r g ks — d s g kr )X k ] = 0. 

This, and an a similar computation that we leave for the reader, proves the 
first part of the proposition. We prove the second part by direct computation: 
Note the following identity: 

[AB, CD] = [A, C]BD + A[B, C]D + C[A, D]B + CA[B, D\. 

Using this identity we find 

[Xr,X'] = \g ri X i , ! f1X j ] 
= [g r \ X.X, + g ri [Xi, g^jXj + g^[g", X,]X t + g sj g ri [Xi, Xj] 
= -gTW^Xj + g^djig^Xi + g sj g ri [Xi, X 3 ] 
= -g"d t (g s ng 3 iX l + g^djig^guXi + g s >g"[X t , Xj] 
= g ri g sj di( 9jl )Xi - g^g^jigu)^ + g s >g"[X t , Xj] 
= g ri g sj (di( gjl )Xi - d 3 {g u )Xi + [Xi,Xj}) 
= g"g SJ F tJ . 

This completes the proof of the proposition. □ 

We now consider the full form of the acceleration terms X k . We have 
already shown that 

didjX k = digj k + djg ik — d k gij- 

Letting 

Fkij = (l/2)(dig jk + djg ik - d k gij), 
we define G r by the formula 

X r = G r + F rs X* + TrstX^XK 

Proposition 6.3 Let T rs t and G r be defined as above. Then both T rst and 
G r commute with the variables Xi. 
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Proof. Since we know that [X;, digjk] = 0, it follows at once that [Xi, F 
0. It remains to examine the commutator [X;,G>]. We have 

[X h G r ] = [X h X r - F rs X* - T rst XsXt] 

= [x h x r ] - [x h F ra x*\ - [x h r rst X'Xt] 

= [X u X r ] - F rs [X h X°] - T rst [X h X»X*\ 
(since F rs and r rst commute with Xi). Note that 

and that 

[Xi,X'X*] = [Xi,X']X* + X'[X h X*] 
= SfX* + X s si. 

Thus 

[X h G r ] = [X h X r ] - F rs 5f - r rst (<JfX* + X'Sft 

= [Xi, X r ] - F r i - YritX 1 - T rs [X s . 

It is easy to see that T^X 1 = T rs [X s . Hence 

[X h G r ] = [X h X r ] — F H — 2T rU X t . 

On the other hand, 

[Xu X r ] = g lr . 

Hence 

[Xi, X r ] = gi r — [Xi, X r \. 

Therefore 

[Xi, G r ] = gi r — [Xi, X r ] — F r i — 2T r ltX t 

= gir - (d r g k i - dig kr )X k - 2r wt X*. 
(since F H = [X r ,X{\ + (d r g k i - dig kr )X k ) Hence 

[X h G r ] = glr - (d r g u - dig*)X % - (dig tr + d t g lr - d^X* 

= gir - (dtgirjX* = 0. 
This completes the proof of the proposition. 
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We now know that G r , F rs and T rst commute with the variables X^. As 
we now shall see, the formula 

X r = G r + F rs X s + T rst X s X t 

allows us to extract these functions from X r by differentiating with respect 
to the dual variables. We already know that 

didjX k = 2T ki j, 

and now note that 

d t (X r ) = [Xi, X r \ = [X t , G r + F rs X° + T rst XsXt] 
= F rs [X t ,X s }+T rst [X l ,X s X i ] 

— Fri ~\~ 21T ritX^ . 

We see now that the decomposition 

X r = G r + F rs X s + TrstX-X' 
of the acceleration is uniquely determined by these conditions. Since 

F rs = [X r ,X s ] = g ri g sj Fij, 

we can regard this result as a description of the motion of the non-commutative 
particle influenced by a scalar field G r , a qauge field F rs , and geodesic mo- 
tion with respect to the Levi-Civita connection corresponding to gij. The 
structural appearance of all of these physical aspects is a mathematical con- 
sequence of the choice of non-commutative framework. 

Remark. It follows from the Jacobi identity that 

Fj 9ir9js-F 

satisfies the equation 

diFjk + djFki + dkFij = 0, 



37 



identifying F^ as a non-commutative analog of a gauge field. Gi is a non- 
commutative analog of a scalar field. The derivation in this section gener- 
alizes the Feynman-Dyson derivation of non-commutative electromagnetism 
[3] where ^ = Sij. In the next section we will say more about the Feynman- 
Dyson result. The results of this section sharpen considerably an approach 
of Tanimura [28]. In Tanimura's paper, normal ordering techniques are used 
to handle the algebra. In the derivation given above, we have used straight 
non-commutative algebra, just as in the original Feynman-Dyson derivation. 

Remark. It is interesting to note that we can rewrite the equation 

X r = G r + F rs X s + r rst x s x* 

as 

X r = G r + [X r , X S ]X* + TsrtX-XK 

(Just substitute the expression for F rs and recollect the terms.) The reader 
may enjoy trying her hand at other ways to reorganize this data. It is im- 
portant to note that in the first form of the equation, the basic terms G r , 
F rs and T rst commute with the coordinates X^. It is this decomposition into 
parts that commute with the coordinates that guides the structure of this 
formula in the non-commutative context. 

7 An Abstract Version of the Feynman-Dyson 
Derivation 

In this section we assume that specific time-varying coordinate elements 
Xi,X 2 ,X 3 of the algebra A are given. We do not assume any commuta- 
tion relations about X±, X 2 , X 3 . We define the field 

H = X xX. 

The field H is an analog of the magnetic field in electromagnetic theory and 
should not be confused with our earlier notation for the Hamiltonian. 
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In this section we no longer avail ourselves of the principle of commuta- 
tivity that is in back of the original Feynman-Dyson derivation. (See the last 
section.) We do not base the derivation to follow on any particular commu- 
tation relations about the variables Xj, but we do take the definitions of the 
derivations that we use from that previous context. Surprisingly, the result 
is very similar to the one of Feynman and Dyson, as we shall see. 

Here A x B is the non-commutative vector cross product: 

(A x B)k = Y%j =1 €ijkAiBj. 

(We will drop this summation sign for vector cross products from now on.) 
Then 

Hk = CijkXiXj = (l/2)eijk[Xi, Xj]. 
We define the field E by the equation 

X = E + X x H. 

We will see that E and H obey a generalization of the Maxwell Equations, 
and that this generalization describes specific discrete models. The reader 
should note that this means that a significant part of the form of electro- 
magnetism is the consequence of choosing three coordinates of space, and 
the definitions of spatial and temporal derivatives with respect to them. The 
background process that is being described is otherwise aribitrary, and yet 
appears to obey physical laws once these choices are made. 

Remarks on the Derivatives. 

1. Since we do not assume that [Xj, Xj] = Sij, nor do we assume [Xj, Xj] = 
0, it will not follow that E and H commute with the Xj. 

2. We continue to define 

di(F) = [F,Xi], 

and the reader should note that, these spatial derivations are no longer 
flat in the sense of section 4 (nor were they in the original Feynman- 
Dyson derivation). 
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3. We define d t = d/dt by the equation 

d t F = F- T,iXidi(F) — F — EiXi[F,Xi] 

for all elements or vectors of elements F. We take this equation as the 
global definition of the temporal partial derivative, even for elements 
that are not commuting with the Xj. This notion of temporal partial 
derivative d t is a least relation that we can write to describe the tem- 
poral relationship of an arbitrary non-commutative vector F and the 
non-commutative coordinate vector X. 

4. In defining 

d t F = F-Z i X i d i (F), 

we are using the definition itself to obtain a notion of the variation of F 
with respect to time. The definition itself creates a distinction between 
space and time in the non-commutative world. 

5. The reader will have no difficulty verifying the following formula: 

d t (FG) = d t (F)G + Fd t (G) + T,idi(F)di(G). 

This formula shows that d t does not satisfy the Leibniz rule in our 
non-commutative context. This is true for the original Feynman-Dyson 
context, and for our generalization of it. All derivations in this theory 
that are defined directly as commutators do satisfy the Leibniz rule. 
Thus d t is an operator in our theory that does not have a representation 
as a commutator. 



6. We define divergence and curl by the equations 
and 

(V x E) k = €i jk di{Ej). 



We now prove a few useful formulas about the vector products. First we 
have the basic identity about the epsilon. 
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Lemma 7.1 Let €ij k be the epsilon tensor taking values 0, 1 and —1 as 
follows: When ijk is a permuation of 123, then e^ k is equal to the sign of 
the permutation. When ijk contains a repetition from {1,2,3}, then the 
value of epsilon is zero. Then e satisfies the following identity in terms of the 
Kronecker delta. 

e abi e cdi = —S ad 5 bc + S ac 5 bd . 

The proof of this identity is left to the reader. The identity itself will be 
referred to as the epsilon identity. The epsilon identity is a key structure in 
the work of this section, and indeed in all formulas involving the vector cross 
product. 

Lemma 7.2 Let A, B, C be vectors of elements of the algebra A. Then 

(Ax B)»C = Ax (B »C). 

Proof. Note that 

(A x B) • C — Tjijj, tij k AiBjC k 

^-'ijk A^tij k BjC k 
= ^ijk AitjkiBjCk 

= A»(BxC). 

This completes the proof of the Lemma. □ 

Lemma 7.3 Let A be any vector of three elements of the algebra A. Then 

VxA = -AxX-XxA. 

Proof. We shall use the summation convention for repeated indices in this 
calculation. 

(V x A) k = eijkdiAj = € ijk [A j: Xi\ 

^ijkAjXi tij k X^Aj 
tjikAjX-i dij^X^Aj 

— -(A x X) k - (X x A) k . 
This completes the proof of the Lemma. □ 
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Lemma 7.4 For A and B any elements in the algebra A, 

V x(Ax B) = -(V • A)B - A* (VB) + (VA) • B + A(V • B) 

where 

(A.(VB)) s = S fc A« 

and 

(((V^.B^E^A)^- 

Proof. 

(V x (Ax B)) c = e abi e cdi d d {A a B b ) 

= (S ad 5 bc + 8 ac 8 w )d d (A a B b ) 
= -S ad S bc d d (A a B b ) + 5 ac 5 bd d d (A a B b ) 
= -d a (A a B c ) + d b (A c B b ) 
= -d a (A a )B c - A a d a (B c ) + d b (A c )B b + A c d b (B b ) 
= [-(V • A)B — Am (VB) + (VA) • B + A(V • B)] c . 
This completes the proof of the Lemma. □ 

Remark. This Lemma, and the observation that the formula in the Lemma 
works in the non-commutative context is due to the author and Keith Bowden 
in conversations around 1999. See [20]. We now give the generalization of 
the Feynman-Dyson result in this formalism. 

Theorem 7.5 With the above definitions of the operators, and taking 
V 2 = dl + d\ + dl H = X x X and E = d t X we have 

1. X = E + X x H 

2. VH = 

3. d t H + V x E = H x H 

4. d t E — V x B = (<9 2 - V 2 )X 
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Remark. Note that this Theorem is a non-trivial generalization of the 
Feynman-Dyson derivation of electromagnetic equations. In the Feynman- 
Dyson case, one assumes that the commutation relations 

[X i ,X j ] = Q 

and 

[Xi,Xj] = 5ij 

are given, and that the principle of commutativity is in place, so that if A and 
B commute with the X^ then A and B commute with each other. One then 
can interpret <9j as a standard derivative with di(Xj) = 5ij. Furthermore, one 
can verify that Ej and Hj both commute with the X { . From this it follows 
that d t (E) and d t (H) have standard intepretations and that H x H = 0. 
The above formulation of the Theorem adds the description of E as dt(X), 
a non-standard use of dt in the original context of Feyman-Dyson, where dt 
would only be defined for those A that commute with JQ. In the same vein, 
the last formula d t E — V x H = (<9 t 2 — V 2 )X gives a way to express the 
remaining Maxwell Equation in the Feynman-Dyson context. 

Proof of Theorem. We begin by calculating 

W = X x H = X x (X x X). 

Hence 

Wi = (X x (X x X))i = -X k (X k X t ) + (X k X t )X k . 
This follows from Lemma 7.1. Hence 

Wi = [X k Xi,X k ] = T, k d k (X k Xi) 

= E fc (9 fc (X fe )Xj + X k d k (Xi). 

But 

dk{X k ) = [X k ,X k ] = 0, 

hence 

Wi = T, k X k d k (Xi) = Xi — d t Xi = Xi — Ei. 

Thus 

X x H = X - E. 
This completes the proof of the first part. 
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V • H = Ei[Hi, Xi] = H • X - X • H 

= (X x X) • X - X • (X x X) = 0, 

since it is easy to verify that (AxB)»C = A»(BxC) for the non-commutative 
vector cross product. 

Since 

d t H = H - X •(VH), 

we have 

d t H + VxE = H-X» (Vff ) + V x E 

= (ixl)-l. (Vff) + V x E 
= XxX + XxX-X» (Vff) + V x E 
= (E + XxH)xX + Xx(E + XxH)-X» (Vff) + Vx£ 
= (ExX + XxE) + (XxH)xX + Xx(XxH)-X» (Vff) + V x E 
= -VxE + (XxH)xX + Xx(XxH)-X» (Vff) + V x E 
= (XxH)xX + Xx(XxH)-X» (Vff) 
= -V x (X x H) - X • (Vff). 
Now, using the formula for V x (A x B), we obtain 

d t H + V x E = (V • X)H + X • (Vff) - (VX) •H- X(V • H) - X • (Vff) 

= (V • X)H - (VX) • H - X(V • if). 
Note that V • X = £j[X;, X] = and that V • if = 0. Thus 

d t H + V xE = -(VX) • if. 

Now note that 

(H x H)k = ((X x X) x H)k = £ijk(X x X)iHj = eijke rs iX r X s Hj 

= (—5j S 5kr + Sj r 5k s )X r X s Hj 
— —XkXjHj + XjXkHj 
= - E,- .V,..V ; 7/ ; = (-(VX) . H) k . 

Therefore 

d t H + V x E = H x H. 
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Now consider 

(V x H)i = -d k (X k X t ) + d k (X t X k ) 
= d k [Xi,X k ] = T, k [[Xi, X k ], X k ] 
= H k d k Xi. 

Hence 

V x H = V 2 X. 

The last part of the Theorem follows immediately from this formula. This 
completes the proof. □ 

Remark. Note the role played by the epsilon tensor e^ fe throughout the 
construction of generalized electromagnetism in this section. The epsilon 
tensor is the structure constant for the Lie algebra of the rotation group 
SO (3). If we replace the epsilon tensor by a structure constant fij k for a 
Lie algebra Qof dimension d such that the tensor is invariant under cyclic 
permutation {fij k = fkij), then most of the work in this section will go over 
to that context. We would then have d operator /variables X i: ■ ■ ■ X d and a 
generalized cross product defined on vectors of length d by the equation 

(A x B) k = f ijk A t Bj. 

The Jacobi identity for the Lie algebra Q implies that this cross product will 
satisfy 

A x (B x C) = (A x B) x C + [B x (A] x C) 

where 

{[B X (A] X C) r = fklrfijkAiBkCj. 

This extension of the Jacobi identity holds as well for the case of non- 
commutative cross product defined by the epsilon tensor. The reader will 
enjoy looking back over this section and seeing that we can still carry Theo- 
rem 7.5 up to the following conclusion with E defined by the second equation 
below. We can no longer take E = d t X, as this depends upon the specific 
properties of the epsilon. 

1. Assume H = X x X. 

2. Assume X = E + X x H. 
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3. Then V • H = 0. 

4. Then d t H + V x E = -V x (X x H) - X • (VH). 

It is therefore of interest to explore the structure of generalized non-commutative 
electromagnetism over other Lie algebras (in the above sense). This will be 
the subject of another paper. 

7.1 The Original Feynman - Dyson Derivation and its 
Gauge Theoretic Context 

The original Feynman-Dyson derivation [3, 11, 9, 23] assumes that we have 
three variables {Xi,X 2 ,X 3 } and the commutation relations 

[Xi,Xj] = 

[Xj, Xj] = 5ij. 

It is also assumed that if A and B commute with the X,, then A and B 
commute with each other. That is, A and B are then "functions of the Xj" . 
We have called this the principle of commutativity. 

With these assumptions one proves that with 

H = X xX, 

(non-commutative vector cross product) and E defined by 

X = E + X x H, 
then E and H satisfy the Maxwell equations in the sense that 

1. E and H commute with the Xj. 

2. V •# = (). 

3. dtH + V x E = 0. 
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where these differential operators have been described in detail (in the non- 
commutative framework) in this section. A key to the original demonstration 
is the principle of commutativity, providing a compass for comparing the re- 
sults with the context of classical calculus. In this section we have seen 
that an abstraction of the Feynman-Dyson argumemt provides a serious gen- 
eralization that encompasses a number of discrete models (to be discussed 
below). In this sub-section, we compare the Feyman-Dyson framework with 
our already-constructed formality of non-commutative gauge theory. 

We use the dynamics 

dX t /dt = X i = P l -A t , 

as before. We restrict to the case where [Xi, Aj] = so that 

9ij — [Xi, Xj] = [Xi, Pj — Aj] = Sij — [Xi, Aj] = 5ij. 

This is the domain to which the original Feynman-Dyson derivation applies. 
We then have 

[Xi,Xj]=0 
[Xi,Xj] = Sij 

and 

R l3 = [Xi, Xj] = diA 3 - djAi + [Ai, Aj\. 

Note that even under these restrictions we are still looking at the possibility 
of a non-abelian gauge field. The pure electromagnetic case is when the 
commutator of Ai and Aj vanishes. In the Feynman-Dyson context, this 
commutator does vanish, since it is given that [Xi, Aj] = for all % and j, 
and the principle of commutativity applies. 

With this interpretation, E is defined by the Lorentz force law 

X = E + X x H 

where H represents the magnetic field. To see how this works, suppose that 
Xi = Ei + FijXj and suppose that E { and F { j commute with X k . Then we 
can compute 

[Xi, Xj] = [Xi, Ej + Fj k X k ] 
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— Fj k [Xi,X k ] — Fj k Sik — Fji. 

This implies that 

F^ = [Xi, Xj] = = diAj — djAi 

since [Xi,XA + [Xj,X,-] = D[X i: Xj] = 0. It is then easy to verify that the 
Lorentz force equation is satisfied with H k = tijkRij and that in this case 
of [Ai, Aj] = leads directly to standard electromagnetic theory when the 
bracket is a Poisson bracket. When this bracket is not zero but the potentials 
Ai are functions only of the Xj we can look at a generalization of gauge theory 
where the non-commutativity comes from internal Lie algebra parameters. 
This shows how a shift of the original Feynman-Dyson derivation supports 
generalizations of classical electromagnetism. 

7.2 Discrete Thoughts 

In the hypotheses of the above Theorem, we are free to take any non- 
commutative world, and the Theorem will satisfied in that world. For ex- 
ample, we can take each Xi to be an arbitary time series of real or complex 
numbers, or bitstrings of zeroes and ones. The global time derivative is 
defined by 

F = J(F'-F) = [F,J], 

where F J = JF' . This is the non-commutative discrete context discussed in 
sections 2 and 3. We will write 

F = JA(F) 

where A(F) denotes the classical discrete derivative 

A(F) — F' — F. 

With this interpretation X is a vector with three real or complex coordinates 
at each time, and 

H = X x X = J 2 A(X') x ApT) 

while 

E = X-Xx(XxX) = J 2 A 2 (X) - J 3 A(X") x (A(X') x A(X)). 
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Note how the non-commutative vector cross products are composed through 
time shifts in this context of temporal sequences of scalars. The advantage of 
the generalization now becomes apparent. We can create very simple models 
of generalized electromagnetism with only the simplest of discrete materials. 
In the case of the model in terms of triples of time series, the generalized 
electromagnetic theory is a theory of measurements of the time series whose 
key quantities are 

A(X') x A(X) 

and 

A(X") x (A(X') x A(X)). 

It is worth noting the forms of the basic derivations in this model. We 
have, assuming that F is a commuting scalar (or vector of scalars) and taking 
A, = X[ - X t , 

di(F) = [F, Xi] = [F, JA t ] = FJA t - JA t F = J(F'A, - A t F) = FA t 
and for the temporal derivative we have 

d t F = J[l - J A' • A]A(F) 

where A = (A 1; A 2 , A 3 ). 

7.2.1 Discrete Classical Electromagnetism 

It is of interest to compare these results with a direct discretization of classical 
electromagnetism. Suppose that X, X', X", X'", • • • is a time series of vectors 
in R 3 (where R denotes the real numbers). Let VX = X' — X be the usual 
discrete derivative (with time step equal to one for convenience). Let A» B 
denote the usual inner product of vectors in three dimesions. 
Assume that there are fields E and H such that 

V 2 X = E + VX x H 

(the Lorentz force law). Assume that E and H are perpendicular to the 
velocity vector VX , and that E is perpendicular to H. 
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Then we have 



VX' x VX = (VX' - VX) x VX = (V 2 X) x (VX) 

= E x VX + (VX x H) x VX 

= ExVX- VX(H • VX) + (VX • VX)H. 

Since E is perpendicular to VX we know there is a A such that E x VX = XH 
and we have H • VX = since H is perpendicular to VX. Therefore 

VX' x VX = XH + | | 2 # 

so that 

H = VX' x £>X/(A+ ||£>X|| 2 ). 
Now using A = A(X), and get 

E = (A' - A) - [A x (A' x A)]/(A + A • A), 

and 

H = (A' x A)/(A + A • A). 

The formula for H is in exactly the same pattern as the formula for H 
in the discrete model for generalized electromagnetism as described in this 
subsection. Up to the time-shifting algebra and a proportionality constant, 
the expressions are the same. The expression for E is similar, but involves 
different time-shift structure. Clearly more work is needed in comparing 
classical discrete electromagnetism with the results of a discrete analysis of 
this generalized Feynman-Dyson derivation. 

7.3 More Discrete Thoughts 

In the Feynman-Dyson derivation of electromagnetic formalism from com- 
mutation relations [11] one uses the relations 

[X t ,X 3 } = 

[Xi,Xj] = kSij 
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where A; is a scalar. In this subsection we shall use 

[Xi,X;]=0 
[Xi,Xj] = JkSij 

as we did in analyzing the one-dimensional case. We shall take 

F = J(F' — F) — [F, J] 

with 

JF = F'J, 

taking the time-step equal to one for convenience. This allows us to have 
scalar evolution of the time series, but changes the issues in the original 
Feynman-Dyson derivation due to presence of the non-commutative opera- 
tor J in the second equation. These issues are handled by the more general 
formalism that we discussed in this section. We aim to see to what extent 
one can make simple models for this version of the Feynman-Dyson rela- 
tions. Models of this sort will be another level of approximation to discrete 
electromagnetism. 

Writing out the commutation relation [X, X] = Jk, and not making any 
assumption that X' commutes with X, we find 

J _1 [X,X] = X'(X' - X) - (X' - X)X 

= X'(X' - X) - X(X' - X) + X(X' - X) - (X' - X)X 
= (X' - Xf + (XX 1 - X'X) = (X' - Xf + [X, X']. 
Thus the commutation relation [X, X] = Jk becomes the equation 

(X 1 -Xf + [X,X'\ = k. 

By a similar calculation, the equation [X, Y] — becomes the equation 

(X' - X)(Y' -Y) + [X, Y'} = 0. 

These equations are impossible to satisfy simultaneously for k ^ if we 
assume that X and X' commute and that X and V commute and that 
[Y, Y] = Jk. For then we would need to solve: 

(X' - Xf = k. 
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(r - Yf = k. 

(X' - X)(Y' -Y) = 0. 

with the first two equations implying that (X — X') and (Y — Y') are each 
non-zero, and the third implying that their product is equal to zero. In other 
words, the equations below cannot be satisfied if the time series are composed 
of commuting scalars. 

[X, X] = Jk 
[Y, Y] = Jk 
[X,Y]=0 

In order to make such models we shall have to introduce non-commutativity 
into the time series themselves. 

Here is an example of such a model. 
Return to the equations 

(X'-X) 2 + [X,X'\ = k. 

(X' - X)(Y' -Y) + [X, Y'} = 

expressing the behaviour for two distinct variables X and Y. If [X, X'] = 0, 
then we have (X' — X) 2 = k so that 

X' = X ± Vk. 

In order for the second equation to be satisfied, we need that 

[X, Y'\ = ±k 

where the ambiguity of sign is linked with the varying signs in the temporal 
behaviour of X and Y. We will make the sign more precise in a moment, but 
the radical part of this suggestion is that for two distinct spatial variables X 
and Y , there will be a commutation relation between one and a time shift of 
the other. 
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If the space variables are labeled Xi, then we can write 



X* +1 = X\ + e\k 

where e™ is plus one or minus one. Thus each space variable performs a walk 
with the fixed step-length k. We shall write informally 

X[ = X, + e,k 

where it is understood that the epsilon without the superscript connotes the 
sign change that occurs in this juncture of the process. We then demand the 
commutation relations 

[X' i ,X j ] = [X' j ,X i ]=e i e j k. 

Each Xi is a scalar in its own domain, but does not commute with the time 
shifts of the other directions. We then can have the full set of commutation 
relations: 

[X' i ,X j ] = [X' j ,X i ]=e i e j k. 
[X i ,X j ]=0 
[X i ,X j ] = Jk8 ij . 

In this system,the elements of a given time series Xi, X[, X", ■ ■ • commute 
with one another. The basic field element in the generalized Feynman-Dyson 
set up is the magnetic field H defined by the (non-commutative) vector cross 
product 

H = (l/k)X x X. 

Here we have 

X t = J{X[ - X,) = JeiVk. 

Thus 

H = J 2 e x e 

where e = (£1,62,63) and e' denotes this vector of signs at the next time 
step. In this way we see that we can think of each spatial coordinate as 
providing a long temporal bit string and the three coordinates together give 
the field in terms of the vector cross product of their temporal cross sections 
at neighboring instants. It is interesting to compare this model with the color 
algebra in [29]. 
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8 The Jacobi Identity and Poisson Brackets 



It is worth thinking through the message of the non-commutative world in 
respect to the existence of the Poisson brackets and their connection with 
continous differentiation and the commutative world of topology and differ- 
ential geometry from which the classical and the quantum models of physics 
are derived. In the classical world there are specific point locations, and the 
notion of a trajectory is given in terms of a continuous sequence of such lo- 
cations. But there is no inherent operational structure intrinsic to the space. 
There is great freedom in the world of commutative and continuous calculus, 
a freedom that allows the construction of many models of temporal evolution. 
Yet we have seen that non- commutative worlds have built in laws, and built 
in patterns of evolution. These patterns of evolution do not lead directly to 
trajectories but rather to patterns of concatenations of operators. At first 
sight it would seem that there could be no real connection between these 
worlds. The Poisson bracket and the reformulation of mechanics in Hamil- 
tonian form shows that this is not so. There is a special non-commutativity 
inherent in the continuous calculus, via the Poisson Bracket. 

It is easy to see the truth of the Jacobi identity for commutators. It is 
just a little harder to see the Jacobi identity of Poisson brackets. It is the 
purpose of this section to recall these verifications and to discuss the nature 
of the identity. 



First let [A, B] = AB - BA. Then 
[[A, B],C] = {AB - BA)C - C{AB - BA) = ABC - BAC - CAB + CBA. 



[[A,B],C] = ABC- BAC 



CAB + CBA, 



[[C,A],B] = CAB- AC B 



BCA + BAC, 



[[B,C],A] = BCA- CBA 



ABC + ACB. 



So 



[[A,B],C] + [[C,A],B] + [[B,C],A]=0. 



This is the Jacobi identity. 
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More generally, a Lie algebra is an algebra A with a (non-associative) 
product ab, not necessarily a commutator, that satisfies 

1. Jacobi identity (ab)c + (bc)a + (ca)b = and 

2. ba = -ab. 

It follows that if we define p a : A — ► A by the equation p a (x) = ax for each 
a in A, then 



so that products go to commutators naturally in the left-regular representa- 
tion of the algebra upon itself. 

Here is another point of view. We have the following equivalent form of the 
Jacobi identity (when ab = —ba for all a and b): 



for all a, x and y in the algebra. This identity says that each element a in 
the algebra acts, by left multiplication, as a derivation on the algebra. In 
this way, we see that Lie algebras are the natural candidates as contexts for 
non- commutative worlds that contain an image of the calculus. 

8.1 Proving the Jacobi Identity for Poisson Brackets 

There are examples of Lie algebras where the non-associative product is not 
a commutator, the most prominent being the Poisson bracket. Here we start 
with a commutative algeb ra C A with two (or more) derivations on CA. Let 



there be operators |~a~ and b | acting on CA (ab is the commutative multipli- 

cation) such that these operators satisfy the Leibniz rule and commute with 
one another: 



Pab = [Pa,Pb], 



a(xy) = (ax)y + x(ay) 





and 




and 



[a 
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for all elements of CA. Then we define the Poisson Bracket on CA by the 
formula 

{a, b}cA = - [b" - ^. 



We wish to see that this product satisfies the Jacobi identity. In order to 
do this we first prove a lemma about the Jacobi identity for commutators 
in a non-associative algebra. We then apply that lemma to the specific non- 
associative product 



Suppose that * denotes a non-commutative and non-associative binary 
operation. We want to determine when the commutator [A, B] = A*B—B*A 
satisfies the Jacobi identity. We first prove a lemma about the Jacobi identity 
for commutators in a non-associative algebra. Let MA be a non-associative 
linear algebra with multiplication denoted by * as above. Let 



and call this the Jacobi sum of a, b and c. We say that the Jacobi identity is 
satisfied for all elements a,b,c G MA if J(a, b, c) =0 for all a, b, c G MA. We 
define the associator of elements a, b, c by the formula 



Let a be an element of the permutaion group S3 on three letters, acting on 
the set {a, b, c}. Let a a , b a , c a be the images of a, b, c under this permutation. 
Let sgn(a) denote the sign of the permutation. 

Lemma 8.1 Let MA be a non-associative algebra as above, then the the 
Jacobi sum J(a, b, c) = [[a, 6], c] + [[c, a], b] + [[b, c], a], for any elements a, b, c G 
A is given by the formula 



Thus the Jacobi identity is satisfied in MA iff the following identity is true 
for all a, b, c G MA. 




J(a,b,c) = [[a,b],c] + [[c,a],b] + [[b,c],a], 



< a,b,c >= (a * b) * c — a * {b * c). 



J(a,b,c) = T lae s 3 sgn(a) < a° ,b° \c a > . 



Z aeS3 sgn(a) <a a ,b a ,c CT >= 0. 
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Proof. For the duration of this proof we shall write ab for a * b. Then 

[[a, b], c] = (ab — ba)c — c(ab — ba) = (ab)c — (ba)c — c(ab) + c(ba), 

[[c, a], b] = (ca — ac)b — b(ca — ac) = (ca)b — (ac)b — b(ca) + b(ac), 
[[b, c], a] = (be — cb)a — a(bc — cb) = (bc)a — (cb)a — a(bc) + a(cb). 



=< a,b,c> — < b,a,c > + < c,a,b > 
— < a,c,b > + < b,c,a > — < c,b,a > . 
This completes the proof. 

Remark. We discovered this lemma in the course of the research for this 
paper. Gregory Wene points out to us that a version of the lemma can be 
found in [22]. We now apply this result to prove that Poisson Brackets satisfy 
the Jacobi identity. 



Theorem 8.2 Let there be operators [a - and b | acting on a commutative 

algebra CA (ab is the commutative multiplication) such that these operators 
satisfy the Leibniz rule and commute with one another: 



Hence 



[[a,b],c] + [[c,a],b] + [[b,c],a] 
(ab)c — (ba)c + (ca)b — (ac)b + (bc)a — (cb)a 
c(ab) + c(ba) — b(ca) + b(ac) — a(bc) + a(cb) 



((ab)c - a(bc)) - ((ba)c - b(ac)) + ((ca)b - c(ab)) 
((ac)b — a(cb)) + ((be) a — b(ca)) — ((cb)a — c(ba)) 





and 




and 



[a 
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for all elements of CA. Define a non-associative algebra MA with product 

a*b = [a - b |. 

Then the commutator in this algebra [a, 6]^ = a*b—b*a will satisfy the Jacobi 
identity. Note that this commutator is the Poisson bracket with respect to 
the above derivations in the original commutative algebra: 

{a, b}cA — [a - b | — [b~ — a] = a*b — b* a = [a, b]j^A- 

This result implies that Poisson brackets satisfy the Jacobi identity. 

Proof. Consider the associator in the non-associative algebra defined in the 
statement of the Theorem: 



< a, 6, c >= (a *b) * c — a * (b * c) = \~&~b\~ c~| — [a" \h~~ c] 



[a - b| — c"| + fa~ b | — c~| — fa" [b~ — c] + \~&\b~ ~~ c] 



= fa - b| — c] + fa~|b ~~ c] 
Note that an expression of the form 



fa" b | — c] 

will return zero when averaged in the summation 

^aes 3 sgn(a) < a a ,b a ,c a > 



since 



fa" b \ ~c] = fa c] b | (the underlying algebra is commutative) and 
these terms will appear with opposite signs in the summation. Therefore we 
find that Jac(a, b, c) = for all a, b, c in R. This completes the proof. 



9 Diagrammatics and the Jacobi Identity 

We have seen that a commutative world equipped with distinct derivations 
that commute with one another is sufficient to produce a non- commutative 
world (via the Poisson brackets) that is strong enough to support our story 
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of physical patterns. Many combinatorial patterns mimic the Jacobi identity, 
and hence provide fuel for further study. In order to illustrate these connec- 
tions, we give in this section a diagrammatic version of the Jacobi identity 
and an interpretation in terms of graph coloring. We will initially work in an 
Lie algebra Q whose product ab satisfies ba = —ab and the Jacobi identity 
a(bc) = (ab)c + a(bc). In Figure 3 we show a diagrammatic interpretation of 
multiplication, consisting in a trivalent vertex labeled with a, b, and ab. As 
one moves around the vertex in the plane, clockwise, one encounters first a, 
then b, and then ab. 




ab = —ba 



ab —ab 

Figure 3 Diagrammatic Multiplication 

In Figure 4 we illustrate the Jacobi identity in the form 

(ac)b = (ab)c — a(bc). 

b c 

v 

q | ^_ 

D" 

Figure 4 — Diagrammatic Jacobi Identity 

To illustrate how this pattern can occur in a different context, consider 
diagrams D of intersecting chords on a circle as shown in Figure 5. By a circle 
we mean a curve in the plane without self-intersections that is a topological 
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circle. By a chord, we mean an arc without self-intersections that is embedded 
in the interior of the circle, touching the circle in two distinct points. Let us 
suppose that we wish to color the chords from a set of q colors such that if 
two chords intersect in an odd number of points, then they receive different 
colors. Let C(D,q) denote the number of distinct colorings of the chords 
of the diagram D, as a function of q. Call such a diagram of intersecting 
chords an intersection graph. We extend such diagrams by allowing internal 
trivalent vertices as illustrated in the abstract by diagram D" in Figure 4 
and by the diagram with the same label, D", in Figure 5. Interpret the 
trivalent vertex as an instruction that all chord lines meeting at a trivalent 
vertex receive the same color. The diagrammatic Jacobi identity of Figure 
4 corresponds directly to the logical coloring identity that says that if we 
have three diagrams D, D', D" with two chords touching in an odd number 
of points in D, one point removed in D', and the two chords fused by a 
trivalent vertex in D" so that they must receive the same color, then the 
number of colorings of D is the number of colorings of D' minus the number 
of colorings of D" . This is just the coloring version of the logical identity 

Different = Anything — Same. 

For graph coloring problems, this identity was first articulated by Hassler 
Whitney [31]. In formulas, we have 

C(D,q)=C(D',q)-C(D",q). 




D D' D" 



Figure 5 — Intersection Graphs 

The convention that we have adopted here - that two chords are colored 
differently if and only if they intersect in an odd number of points, makes a 
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demand on the interpretation of the trivalent nodes. All arcs entering a given 
node must receive the same color. After more nodes are added we will have 
connected components of the resulting graph that contain nodes (the outer 
circle is not regarded as part of the graph). Call such a connected component 
a web in a given diagram. Each web is colored by a single color. We regard a 
chord without nodes as a (degenerate) web. We take the convention that if 
the total number of intersections between two distinct webs is odd, then they 
must receive different colors. Of course, a web may have self-intersections; we 
define the sign of the coloring of a given web to be —1 if it has an odd number 
of self-intersections and +1 if it has an even number of self-intersections. The 
sign of the coloring of a diagram is the product of the signs of its component 
webs. Note the the sign of a chord is positive. With these conventions, 
the formulas in Figures 4 and 5 match perfectly and can be understood as 
indicating parts of larger diagrams that differ only as indicated. We see, as in 
Figure 6, that an extra self-intersection added to a trivalent vertex changes 
the sign of its web. This corresponds to the algebraic interpretation of such 
as vertex as ab = —ba. See Figure 3. 
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Figure 6 — Verifying the Twist Identity for Color Diagrams 

In Figure 6 we illustrate how these sign conventions are consistent with 
the coloring formula/ Jacobi identity. In this figure, we begin with the Jacobi 
identity with a twist (crossing) added to each diagram. The original diagram 
with one crossing now has two, and hence is equivalent to a diagram with none 
(no local requirement of difference). The original diagram with no crossing 
now has one, and is interpreted as a requirement of difference. Rearranging, 
we find the Jacobi identity again, but with an extra crossing and change of 
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sign for the noded diagram. The conclusion is that adding a crossing to a 
node changes the sign of its diagram. 



We see that the patterns of counting colorings of chord diagrams corre- 
spond formally to the axioms for a Lie algebra. This example indicates how 
a combinatorial context can lead to the very formalism on which this paper 
is based, but though different structures than one could have initially visual- 
ized. Diagrammatic Lie algebras similar to this example feature prominently 
in the theory of Vassiliev invariants [1, 21] of knots and links, and may form 
the basis for new models for the structures that we have discussed in this 
paper. 



We have concentrated in this section on a coloring example not only 
because the occurrence of the Jacobi identity in this context may appear 
startling, but also because there is a more direct relationship with coloring 
in regard to the fundamental Lie algebra of SU(2) (or equivalently SO (3)) 
that underlies the structures we have discussed in this paper. The Lie al- 
gebra of 5*0(3) has structure constant the alternating epsilon symbol 
that we have used again and again in Section 7 for the generalization of the 
Feynman-Dyson derivation. This epsilon can be expressed diagrammatically 
as a trivalent vertex. The basic epsilon identity 



Cabi^cdi — —8ad8bc + &achd 



can be written diagrammatically, and it leads at once to a diagrammatic Ja- 
cobi identity. See Figure 7 for the diagrammatic form of the epsilon identity. 
The epsilon itself is closely related to coloring (See [27, 4, 10]), but that is 
another story and we shall stop here. 



63 



a b 




^abi^cdi $ac$bd $ad$bc 

Figure 7 — The Epsilon Identity in Diagrammatic Form. 

10 Epilogue 

We have sought in this paper, to begin in an algebraic framework that nat- 
urally contains the formalism of the calculus, but not its notions of limits or 
constructions of spaces with specific locations, points and trajectories. It is 
remarkable that so many patterns of physical law fit so well such an abstract 
framework. We believe that this is indicative of the secondary nature of 
point sets, topologies and classical differential geometries in physics (Com- 
pare [2]). In this paper we have dispensed with spacetime and replaced it 
by algebraic structure. But behind that structure, the space stands ready to 
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be constructed, by discrete derivatives and patterns of steps, or by starting 
with a discrete pattern in the form of a diagram, a network, a lattice, a knot, 
or a simplicial complex, and elaborating that structure until the specificity 
of spatio-temporal locations appear. 

There are many ideas for producing location. Poisson brackets allow us to 
connect classical notions of location with the non-commutative algebra used 
herein. Below the level of the Poisson brakets is a treatment of processes and 
operators as though they were variables in the same context as the variables 
in the classical calculus. In different degrees we have let go of the notion of 
classical variables and yet retained their form, as we made a descent into the 
discrete. 

In order for locations to appear from process, one may need an appropri- 
ate degree of recursiveness. Lie algebras begin the process with their fully 
self-operant structure of derivations. It is just such bootstrapping that fits 
into the basis of our concerns and produces the ways to make spaces emerge, 
through process, from abstract algebra. 
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